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1. Introduction 

Various nonlinear dispersive wave systems exhibit instability, known as the modulational 
instability (MI). MI stems from the interaction between nonlinearity and group- velocity 
dispersion and manifests itself as a self-induced amphtude modulation of a continuous 
wave propagating in a nonlinear medium, with subsequent generation of localized 
structures. MI arises in many physical settings including plasma [H [21 [3l H] , fluids [5], 
nonlinear optics [HI [3 El E] , and atomic Bose-Einstein condensates [TOl ITT] . 

In the context of nonlinear fibre optics, MI requires anomalous dispersion of a 
medium to exist. In this case, negative group- velocity dispersion combined with the 
self-phase modulation amplifies the modulational frequency sidebands producing a train 
of ultra-short pulses as a result of breakup of a continuous wave [121 El [13] . 

On the other hand, when two or more optical beams propagate down the fibre, 
MI can occur even in the normal- dispersion regime, at the cost of the cross-phase 
modulation [TH [15]. Such a situation is usually described by a set of coupled nonlinear 
Schrodinger (NLS) equations. A thorough analysis of MI for two nonlinearly interacting 
waves in a Kerr medium was performed in recent papers [T6l [T7] . 

The above results concern the class of local nonlinear media when the response 
of the medium at a particular point depends solely on the wave intensity at 
that point. Alternatively, nonlocality plays a key role in physical systems where 
transport phenomena and finite-range interaction cannot be neglected. The examples 
of importance of spatial nonlocality for the development of MI can be found in 
plasmas [TH [19], hquid crystals [20], discrete nonlinear systems [21], Bose-Einstein 
condensates [22]. In particular, nonlocality changes drastically the growth rate and 
bandwidth of instability caused by stochasticity of parameters of a Kerr medium |23j . 
The present status of MI of scalar waves in nonlocal media is discussed in the review 
papers [MI25]. 

In our paper we extend the analysis of MI in coupled wave systems to the case 
of nonlocal media. It was demonstrated recently that vector soliton structures in 
nonlocal media exhibit properties that have no counterpart in the scalar case |26j . 
In particular, a possibility was revealed to stabilize multipole (dipole, quadrupole, 
hexapole, etc.) vector solitons in bulk nonlocal media at the cost of vectorial coupling 
|27] . In this case mutually incoherent nodeless and multipole components jointly 
induce a nonlinear refractive index profile. As regards higher-order nonlinearities, a 
one-dimensional phenomenological model of a nonlocal medium with focusing cubic 
and defocusing quintic nonlinearities was developed in [28]. Two types of solitons, 
even-parity (fundamental) and odd-parity (dipole) were found and their stability was 
explored. Note that solitons of this sort are unstable in local media. Lastly, an 
interesting result should be mentioned concerning formation of stable three-dimensional 
spatiotemporal solitons in nonlocal Kerr media [29]. Fundamental soliton of this model 
is stable provided its propagation constant exceeds a certain critical value inverse 
proportional to a nonlocality parameter. 
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Below we derive the dispersion relation which characterizes the propagation 
of two interacting waves in a nonlocal Kerr-type medium. For the case of zero 
group-velocity mismatch we are able to perform a complete analytic description of 
stability/instability regimes for generic wave amplitudes. It is shown that nonlocality 
suppresses considerably the growth rate and bandwidth of instability. For nonzero 
velocity mismatch, we investigate a role of nonlocality in a management of the 
stability properties in the system under consideration. Nonlocality is shown to provide 
an efficient stabilization mechanism permitting stabilize interacting waves which are 
strongly unstable in a local medium. 



2. Model 



Propagation of two nonlinearly interacting waves with complex amplitudes u{t,x) and 
v{t, x) in a nonlocal Kerr medium is modelled by the following system of dimensionless 
equations: 

f du ,,du\ did'^u 



^dt dx J 2 dx'^ 

+siu J dx'R{x-x') (\u{t,x')\'^ + a\v{t,x')\' 

' dv ,,dv\ dod'^v 
+ V2— 1 + 



(1) 



\dt dx J 2 dx"^ 

+S2V J dx'R{x-x') {\v{t,x')\'^ + a\u{t,x')\'^^ =0, 

which naturally generalize those [SD] for a local medium. 

Physical meaning of coordinates t and x depends on the context of the problem: 
they are time and coordinate for Bose-Einstein condensates, longitudinal coordinate 
z and transversal coordinate x for optical beams in waveguides, the propagation 
coordinate z and time t for optical pulses in fibres, respectively. The parameters dj 
(j = 1, 2) determine effective atom masses in condensate models, diffraction coefficients 
in waveguides and group-velocity dispersion for pulses. The parameters Vj stand for 
group velocities, Sj determine nonlinearity coefficients, and a > is a cross-phase 
modulation parameter. No restrictions are posed on signs and magnitudes of dj and Sj. 
R{x) is a positively definite symmetric response function of a nonlinear medium obeying 
the normalization condition / dxR{x) = 1 (we do not specify the integration limits 
when the integration is carried out along the whole line). An important development of 
nonlocal models with asymmetric response functions relevant, e.g., to account for the 
Raman effect, was achieved by Wyller [31]. 

We will illustrate our results by the exponential response function 

Here a is a nonlocality parameter. The function ([2]) occurs as a solution of the diffusion- 
type equation for nonlinearity [26] and adequately describes the nonlinear response of 




Modulational instability and nonlocality management in coupled NLS system 



4 



thermo-optical materials, liquid crystals and partially ionized plasmas. For cr — we 
have R{x) 5{x) and hence reproduce the local model used in |17j . 



3. Modulational instability: general analysis 

In order to investigate the modulational instability in the system ([1]), we start with the 
continuous wave solution of the form 

u = Moe"'^^*, V = t;oe-'"2*, (3) 

where Uq and Vq are constant real amplitudes, and frequencies uj are given by the 
relations 

ui = -si (ul + avl'j , UJ2 = -82 {aul + vf^ , 

which are the same as for the local coupled NLS system. 

Following the standard procedure, let us consider a small perturbation of the 
continuous waves ([3]): 

u={uq + a{t, x))e~'"^'\ v = {vo + b{t, x))e-''^'\ (4) 

Substituting (jll) into ([T]) and neglecting terms with the second and higher orders of the 
complex perturbation amplitudes a and b, we obtain a system of linear equations 
. f da da\ di d^a 



+ V,-£-j + y ^ + siul J dx'R{x - x'){a + a*){t, x') 
+asiUoVo / dx' R{x — x'){b + b*){t, x') = 0, 



dt 



+as2UoVo J dx'R{x — x'){a + a*){t,x') = 0. 

We seek for solutions of ([5]) in the form 

a = A+e^^'^^~^*^ + A* Q-^('^^-^*'t) ^ 

where k and Q are the wave number and frequency of the perturbation, respectively. 
Hence, the condition of solvability of the system ([5]) takes the form of the dispersion 
relation 

\n - ViKf - nj] [{n - V2Kf - nj] - ^^2 = o, (6) 

where 

Ql = ^diK^ {diK^ - 4siMo-R) , VlI = jd2K^ {d2K^ - 4s2fo-R) ) 

^12 = did2SiS2 (auoVoK,"^ Rj 
Note the appearance of the Fourier transform R{k,) of the response function R{x), 

RiK) = [ dxR{x)e''''^ = \ ^ , Ri-n) = R*{k), (8) 



4' V- ' 4- V- - ' (7) 
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in the expressions for Q'j and fifj- Since R{k) — > 1 for a ^ 0, the expression (jS]) reduces 
in this hmit to the well-known dispersion relation for the local model fi5[ [T6| [T7]. 

The dispersion relation ([6]) takes the form of a general fourth-order algebraic 
equation for Q. In order to provide better insight into the properties of its solutions, 
we begin the analysis from the particular case of equal group velocities, Vi = V2, which 
allows us to perform the complete analytical treatment. 



4. Equal group velocities 

If Vi = V2 = V, we can eliminate the terms with Vj from ([5]) by the Galilean 
transformation t' = t, x' = x — Vt, and the dispersion relation ([6]) reduces to the 
biquadratic polynomial in Q, 

- Ql) (q^ - Ql) - = 0. (9) 



Introducing the definitions 

N = Sjdj, ip^ = ^{dl + dl), ^^^^ 

we write equation ^ in the form 

^2 (11) 

— Idfdl - AN (j^idlul + I22dlvl) + 16iV^(l - a'^)i2ii22ulv^ = 0. 
lb L V / J 

In the local models, stability properties are determined directly by the analysis of the 
function Q{k,). In contrast, to establish general properties of the stability /instability 
domains for the nonlocal case, we will primarily consider the function Q{N). More 
detailed information will be then provided by a particular choice of the positively definite 
function N{k). 

Solving ffTTj) gives 



where 



{cp' - iV^ ± v^) , (12) 
D{N) = AAiN"^ + A57]N + 5^ (13) 



S=-(dl-dl 



(14) 



and 

Ai = nifi2{auoVoY + ^7^ ^ ^ ^ (^^^o ~ ^^^0] 
1 
A 

Without loss of generality we take S > 0. 

Stability of continuous waves ([3]) will be ensured if both Q'^ and Q"^ are positive 

for any choice of the continuous wave amplitudes uq and vq. In other words, we do not 

consider situations when stability/instability conditions fulfil only for specific relations 

between uq and Vq, because they are not generally satisfied. Evidently, we get fl^. > 
if D > and \/D < _ jy^ 
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Analysis 

To reveal positiveness regions of we analyze at first the function D{N) ( fT3l) . 
Positiveness of D{N) depends on Ai and 77. It is seen from ( fT4l) that Ai > for 
any a, Uq and Vq if /ii/i2 > 0, as well as for a < 1 if /ii/i2 < 0. Note that Ai could be 
positive for /ii/i2 < and a > 1 under a specific choice of uq and vq, but we disregard 
such a possibility, as was stated above. Hence, in what follows we consider Ai to be 
positive for generic choice of uq and vq, and for only restriction a < 1 if /ii/i2 < 0. 

As regards the role of rj, it can be elucidated from the analysis of the discriminant 
Di = — 16/ii/i2(a5'Uo'yo)^ of the equation D{N) = 0. Evidently, D > ii Di < 0, i.e., 
for fiifi2 > and arbitrary a (the parabola D{N) = lies above the horizontal axis A^). 
If > that occurs for fiifi2 < and hence for a < 1 (two real intersection points of 
the parabola with the axis), we have two possibilities for D > 0, depending on the 
sign of 7]-. 

(i) 1] > means /ii < and /i2 < 0. Both intersection points lie on the negative 
(unphysical) part of the A^ axis, and D > for any A^ > 0; 



(ii) T] <0 means < and fj.2 > 0. Hence, D > for A^ G (0, N'_) U(A'|, 00), 



Lastly, D < if > 0, a < 1, ?7 < 0, and A^ G {N'_,N'^). Let us summarize the 
information obtained from the above analysis for the generic uq and vq'- 
D > for Ai > and 

(i) /ii/i2 > 0, V(a, A^); 

(ii) yUi > 0, /i2 < 0, a < 1, VA^; 

(iii) /xi < 0, /i2 > 0, a < 1, AT G (0, A^:) U(A^;, 00); 

D < for Ai > and 

(iv) /ii < 0, /i2 > 0, a < 1, A^ G {N'_, N'^). 

To study the condition ^/D < ip^ — Nil>, we introduce Z = D — [ip^ — NipY, or 



where 




(15) 



Z = 4A2A^^ + 4/A^ + g, 



(16) 



where 





g = -djdl. 



They are negative values of Z that correspond to the condition \/D < ip^ — Nij), and 
we require ipP' > Nip. 
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The discriminant D2 = 16(/^ — A25') of the equation Z{N) = has the form 
D2 = {jiid\ul + li2dlv'^ - 4(1 - a;^)/ii/i2 {did2UQVQf . 

It is seen that < occurs only for special requirements on uq and Vq] hence, we take 
in what follows -D2 > for generic Uq and Vq. This condition takes place for 111112 > ^ 
and any a, as well as for /ii/i2 < and a < 1. 

Now we can establish properties of Z{N). First of all we see that the presence of 
/ in ( fT6l) excludes the case /ii/i2 < because it would lead to nongeneric conditions on 
uq and vq. Solutions of Z = are written as 

Since D2 > 0, both Z+ and are real. At first we consider the case A2 > which 
occurs for /ii/i2 > and a > 1. If / > 0, i.e. /ii > 0, ^2 > 0, and hence a > 1, we 
obtain > and Z_ < 0. Therefore, Z will be positive for > Z^ and negative for 
G (0, For / < 0, i.e., /ii < 0, /i2 < and a > 1, the solutions of Z{N) = 0, 

= ^ {±^P + \dldlA2 + I/I j (17) 

obey the same properties as Z±, i.e., Z > for A^ > and Z < for A^ G (0, Z'_^). 
If A2 < (i.e., fiifj.2 > and a < 1), both roots of Z{N) = are positive for / > 

{fll > 0, /i2 > 0), 

^± = ^ {dlul + dy, ± v'(« + «'-|A2M?t/i) > 0, 
and negative for / < (/ii < and /i2 < 0), 

Z\ = {dlul + dlvl ± ^ {dlul + dlvlf - \A2\dldij < 0. 

Hence, for A2 < we have Z{N) > for / > and A^ G (Z_, while Z < for 
/ > and A^ G (0, Z_) U(^+j 00), as well as for / < and any A^. 

In addition, we should provide the condition ip"^ > Nip to fulfil. Two possibilities 
exist to satisfy this condition, depending on the signs of /ij: 

(i) for yUi > and jU2 > and any a we should assume the restriction N < B to fulfil, 
where 

1 d"^ ^ d"^ 

B = \ ^^"^^ (18) 

(ii) there are no restrictions on A^ when /ii < and /i2 < for generic and fo and 
any a. 

On the other hand, < A^?/^ for /ii > 0, //2 > and N > B. 
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Figure 1. Maximum wave number Kmax of the instability range (0, Kmax) versus 
the nonlocality parameter a. Nonlocality suppresses the instability bandwidth. The 
instability regime corresponds to the case (g) in ([^U]) . 



4-2. Summary of the analysis 

Summarizing the results obtained from the analysis of solutions Vl\ f|T2|) . we infer that 
stability regions {D > 0, ^/D < (p^ — N^!) are determined by the following conditions 
on the parameters of equations ([T]) for generic choice of the amplitudes Uq and Vq\ 

(a) yUi < 0, /i2 < 0, a < 1, VA^; 

(b) /ii > 0, /i2 > 0, a > 1, e (0, min(5, 

(c) /ii < 0, /i2 < 0, a > 1, AT G (0, Z;); (19) 

> > a < 1 I ^ ^ (0,min(i?,Z_)), B < 
' '[Ne{0,Z^)[j{Z\,B), ifZ\<B. 

Accordingly, instability regions are characterized by the parameters 

(e) fxi <0, fX2> 0, a <l, N E {N'_, A^+); 

(f) yUi > 0, /i2 > 0, a > 1, > min(i?, Z+); 

(g) /ii < 0, ^2 < 0, a > 1, A^ > ^^^^ 

r N e {Z^,Z+), HB <z^, 

(h) /ii > 0, /i2 > 0, a < 1, I N e (B, Z+), if 5 G Z+), 

[ N > B, if B> Z+. 

It should be noted that stability properties depend on the sign of the products 
fij = Sjdj and not on Sj and dj separately. Note also that we reproduce the results of 
Kourakis and Shukla [T7j in the local limit (A^ = /t^^). At the same time, the availability 
of an additional degree of freedom in the form of nonlocality enables us to manage the 
conditions of stability/instability. In the next subsection we demonstrate this possibility. 
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4-3. Examples 

The instability growth rate is characterized by the quantity 



A(k) = max 



im Jni 



im Jni 



(21) 



Below we use the exponential response function ([2]). For this function we can invert 
the dependence N{k): 

We fix magnitudes of all the parameters, except for a, to reveal the role of nonlocality 
in the development of MI for a number of regimes listed in (I20p . As regards the cross- 
phase modulation parameter a, we take its featured values, namely, a = 2 for a > 1 
and a = 2/3 for a < 1. At first we consider the regime (g) in ( l20l) and choose the 
parameters as fii = fi2 = Uq = Vq = 1, df = d^ = 4, a = 2. They correspond 
to the defocusing nonlinearity in planar waveguides or normal dispersion in a focusing 
Kerr medium (recall that MI does not occur in such conditions for a single field). Since 
N > Z', for the regime (g), we have k G (0, Uraax) from f|T7j) with 



"'max 



1 + ^-1 



2a2 

Figure [1] shows the dependence fi;max(c") with Z'^ from f|T7|) . We see that the range 
(0, Kmax) of instability wave numbers narrows with growing nonlocality. Figure [2] shows 
the MI gain spectra A(/t) f l2Tl) for three values of the nonlocality parameter a. The 
maximum gain Amax decreases with an increase in the nonlocality parameter cr. Such a 
decrease is illustrated by figure [3l 

The second example covers the regime (e) in fl20|) . This regime describes a joint 
propagation of two waves, when one of them 'sees' the normal dispersion of a medium, 
while the second wave moves in the anomalous dispersion environment. Since the 
allowable values of lie within the interval {N'_,N'^), we have two limiting values 
of the modulation wave numbers, namely, K- and 



4a2 
1 + ^-1 



2a2 \\\ N!^ 

We choose /ii = —1, /i2 = 1, Uq = Vq = 1, df = 2, d"^ = 4:, a = 2/3. As it is seen 
from figure m the interval width {k,_,k^) decreases with increasing a. Owing to the 
dependence of k± on the nonlocality parameter a, the MI gain spectra demonstrate 
different boundary values of the bandwidth position, both minimal and maximal, for 
different a (figure [5]). 

5. Group-velocity mismatch and nonlocality management 



When we account for group- velocity mismatch (Vi ^ V2), the dispersion relation ([6]) 
gives a general fourth-order algebraic equation. To avoid cumbersome calculation, we 
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Figure 3. Dependence of the MI gain maximum Amax on the nonlocaUty parameter 
a for the instability regime (g). 




Figure 4. Dependence of limiting values k± of modulation wave numbers on the 
nonlocality parameter a for the instability regime (e). The instability bandwidth 
shrinks with increasing nonlocality. 
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Figure 5. MI gain spectra for tlie instability regime (e) for different values of 
nonlocality parameter a. Three curves correspond (from right to left) to cr = 1, 2, 
3. 

will follow [321 [T7] and give a geometrical description of the nonlocality management in 
the case of the group-velocity mismatch. 

Let us write the dispersion relation ([6]) in the form 

(f] _ v^Kf -nl = -. (23) 

This formula describes a mutual arrangement of two curves - a parabola on the left 
and a more complicated 'right' curve. We know from the previous exposition that the 
stability regime is provided by four real solutions of the dispersion equation ([6]), or, in 
other words, by four intersection points of the above curves. If the number of intersection 
points is less than four, the system is unstable. We show below that nonlocality of a 
medium provides an efficient tool to manage stability/instability properties giving a 
possibility to assure stability of a system which is unstable in a local medium. 

Figure [6] corresponds to light wave propagation in a local medium (cr = 0) with the 
parameters /ii = /i2 = 1, uq = = 1, d\ = d"^ = a = 2/3, = 2, and Vi = V2 = 0, 
i.e., zero group- velocity mismatch. The choice of parameters gives fi^ = fi^ = 2 and 
^2^2 = 16/9 and means individual stability of both waves. There are four real intersection 
points of the parabola with the 'right' curve, the latter consists of two hyperbolae and 
a curve in the lower half plane with two singular points. Therefore, this situation 
corresponds to absolute stability. Let us now introduce the velocity mismatch, taking 
Vi = and V2 = V2. The 'right' curve shifts horizontally, losing two intersection points 
and thereby producing instability (figure [7]). If, however, nonlocality is coming into 
play [we take cr = 1 and recalculate Qf, and fifg accordance with ([7])], the four 
intersection points occur again, and stability is restored (figure [H]). 

Our second example is concerned with a stable-unstable wave pair determined by 
the parameters /ii = —1.5, fi2 = ^, Uq = vq = 0.6, a = 2/3, = 1 which correspond 
to f2f = 1.1025, ^2 = —0.11 and ^2^2 = —0.0864. For the group- velocity mismatch 
(Vi = 0.2, V2 = 0) there are only two intersection points - the system is unstable (figure 
|9]). Retaining this mismatch but allowing nonzero nonlocality (a = 0.5), we restore the 
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Figure 6. Zero group-velocity mismatch {Vi — V2): four intersection points of the 
curves (the parabola is depicted by dashed line) defined in (13) provide stability of the 
continuous wave solution for a local medium {a = 0) with the parameters d\ = = 4, 
Ail = = 1, = 2,ul = vl = 1. 
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Figure 7. Breaking of stability (two intersection points only) in a local medium 
(ct — 0) caused by nonzero group- velocity mismatch (Vi = 0, V2 = V2). Other 
parameters are the same as in figure [6l 




Figure 8. Restoration of stability (four intersection points) in a nonlocal medium 
[a — 1) with nonzero group- velocity mismatch (Vi — 0, V2 = \/2). Other parameters 
are the same as in figure [HI 
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Figure 9. Stable-unstable wave pair in a local medium with the group-velocity 
mismatch. The parameters are jii — —1.5, /i2 = 1, ttg = ^ O-^j — 2/3, = 1, 
Vi = 0.2, V2 = 0. The system is unstable because there are only two intersection 
points of both curves. 
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Figure 10. Restoration of stability of the stable-unstable wave pair in a nonlocal 
medium (cr — 0.5). Other parameters are the same as in figure [51 



system stability (figure [TOl) . 

The last example provides the most emphatic illustration of the nonlocality 
management. Let us consider the unstable- unstable wave pair which corresponds to 
negative values of Vl{ and Vt\. It is stated in [17] that such a pair is always unstable 
in a local medium (see figure [TTl no intersection points). The situation is drastically 
altered in a nonlocal medium. It is seen from (I7l) that fi^ can change sign for sufficiently 
large nonlocality parameter cr, providing stability of the system. This possibility is 
exemplified by figure [121 

6. Conclusion 

We have investigated the MI of two interacting waves in a nonlocal Kerr medium, 
both for zero and nonzero group- velocity mismatch. The case of zero mismatch admits 
the full analytical description of stable/unstable regimes of the wave propagation. We 
have demonstrated that nonlocality suppresses the growth rate peak and bandwidth of 
instability. For the case of nonzero velocity mismatch we have revealed a crucial role of 
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Figure 11. Unstable- unstable wave pair in a local medium. The parameters are 
^1 = /i2 = 0.3, uo = vo = 0.6, a = 2/3, = 1, Vi = 0.1, V2 — Q which correspond 
to Vl\ = = —0.0855. There are no intersection points, the system is absolutely 
unstable. 
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Figure 12. Unstable-unstable wave pair in a nonlocal medium {a = 3, Vl\ 
0.0117). Other parameters are the same as in figure [TT] The system is stable at the 
cost of nonlocality management. 



nonlocality in management of the stability properties in the system under consideration. 
In particular, we can provide stability of the nonlocal system in a regime when its local 
counterpart demonstrates absolute instability. 
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